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ABSTRACT 

The Galactic Center lobe is a degree-tall shell seen in radio continuum images of the Galactic center 
(GC) region. If it is actually located in the GC region, formation models would require massive energy 
input (e.g., starburst or jet) to create it. At present, observations have not strongly constrained the 
location or physical conditions of the GC lobe. This paper describes the analysis of new and archival 
single-dish observations of radio recombination lines toward this enigmatic object. The observations 
find that the ionized gas has a morphology similar to the radio continuum emission, suggesting that 
they are associated. We study averages of several transitions from H106a to H191e and find that the 
line ratios are most consistent with gas in local thermodynamic equilibrium. The radio recombination 
line widths are remarkably narrow, constraining the typical electron temperature to be less than 
about 4000 K. These observations also find evidence of pressure broadening in the higher electronic 
states, implying a gas density of n e = 9IOI45Q cm -3 . The electron temperature, gas pressure, and 
morphology are all consistent with the idea that the GC lobe is located in the GC region. If so, the 
ionized gas appears to form a shell surrounding the central 100 parsecs of the galaxy with a mass of 
roughly 10 5 Mq, similar to ionized outflows seen in dwarf starbursts. 
Subject headings: Galaxy: center — radio line: general 



1. INTRODUCTION 

Sinc e their pred i ction a nd discovery more than 4 years 
ago ( Kardashcv 1959; D ravskikh fe Dravskikhl Il964t 
ISorochenko fe Borozichlll964h . radio recombination lines 
have been useful tools for probing the physical condi- 
tions of the ISM. The properties of thermal emission can 
constrain fundamental parameters such as ga s tempera- 
ture and density (e.g., lRoelfsema fc G oss 1992). Further- 
more, the spectral line can provide valuable information 
about the kinematics of the emitting gas. 

In the Galactic Center (GC) region, there have been 
extensive observations of radio recombination lines. 
These lines are useful because they are detectable de- 
spite heavy extinction and can be used to separate 
ther mal H II regions from the widespread nonthermal 
gas ()Pauls fc Mezgerl [1975). Separating thermal and 
nonthermal emission is especially important in the GC 
region, which is filled with extende d ionized gas and 
synchrotron-emitting features (e.g.. lYusef-Zadeh et al.1 
[l98llMorris fc SerabvdlilM ILang et al.lll997l) . Recom- 
bination line surveys have found that diffuse ionized gas 
is strongly correlated with the molecular gas density and 
star formation, suggesti ng that the GC regio n should be 
filled with ionized gas (|Lockman et al.l H996). This gas 
can provide a direct probe of the conditions in the GC 
region. 

The information available by studying ionized gas 
could help address questions surrounding a feature in 
the GC region known as the Galactic center lobe (GCL) . 
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The GCL was first noted in a 10 GHz radio contin- 
uum survey of the ce ntral few degrees of the Galaxy 
(|Sofue fc Handalfl98l . The survey found a degree-tall 
loop of emission north of the galactic plane, surrounding 
the central degree, covering an area from I = 359? 2 — 0?2 
and b = 0?2 - 1?2. The resemblance of the GCL 
to a wind-driven mass outflow was str iking and led to 
several models for it s formation (e.g., [Chevalier 1992; 
iMelia fc Fal"ckell2001l ). However, while the first study of 
the radio spectral index foun d it to be consistent with 
thermal emission (lSofudl!985h . later work suggested the 
GCL was nonthermal dTsuboi et alj 119861 ; iReich et"ai1 
119871) . iBland-Hawthorn fc Cohen! (|2003l ) noted the mor- 
phological connection between the radio continuum and 
mid-IR emission and proposed a starburst outflow model. 
However, the connection of the ionized gas to this feature 
is still unclear. 

Motivated by a desire to understand the nature of the 
emission in the GCL, we have conducted new radio re- 
combination line observations and collected archived ob- 
servations in this region. This paper describes observa- 
tions of the GCL with the Robert C. Byrd 100-m Green 
Bank Telescope 5 (GBT) and the Hat Creek Radio Ob- 
servatory 26-m telescope (HCRO). Section [5] shows how 
the observations were made and calibrated. Section[3]de- 
scribes the results of the observations, including the mor- 
phology, line characteristics, and kinematics of the line 
emission in the region. The HCRO observations show 
that the HI 09a emission has a morphology similar to 
the radio continuum emission. The GBT observations 
find that the recombination line emission is this region 
is in local thermodynamic equilibrium (LTE). The un- 
usually narrow recombination line emission throughout 
the GCL places useful constraints on the gas tempera- 

5 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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ture and other properties. In § |H we discuss the derived 
quantities of the ionized gas, including gas density, tem- 
perature, and the total mass (assuming a geometry and 
distance) . Section [5] summarizes the results of this pa- 
per. A detailed discussion of the nature of the GCL is 
left for a later paper (Law, in prep). 

2. OBSERVATIONS AND DATA REDUCTIONS 

2.1. HCRO 

In 1985, we used the HCRO 26 m telescope to map the 
H109a recombination line toward the GCL. At 5 GHz, 
the HCRO beam size is 10' and its Nyquist sample size 
is 3'.3. Pointings were made over a roughly 1° square 
region on a 6' grid; neighboring points in the grid were 
skipped to form a checkerboard pattern. Spectra were 
calibrated by position switching to a distant, fixed posi- 
tion, with a typical integration time of 2 hours for each 
spectrum. Several scattered points outside the radio con- 
tinuum emission associated with the GCL were also ob- 
served to constrain the background emission. A crude 
gain calibration was made by comparing the brightness 
temperatures between the HCRO and GBT observations 
described below; this comparison implies a gain of ~ 13 
JyK" 1 . 

The spectrometer at HCRO was tuned to a central 
frequency of 5008.923 MHz, 512 channels, and a 20 
MHz bandwith. This bandwidth covered the H109a 
and H137/3 transitions, although only the H109a tran- 
sition was studied. Spectra were Hanning smoothed and 
a seventh-order baseline fit before Gaussian line fitting. 
The velocity uncertainty of the fits are about 1 km s _1 . 
The system temperature was assumed to be 40 K, based 
on observations of empty sky. 

Figure [T] shows images of the H109a line antenna tem- 
perature and velocity from the HCRO data. The line 
brightness is not corrected for atmospheric absorption. 
Since the checkerboard observing pattern makes visual- 
ization awkward, the data were interpolated onto a regu- 
lar grid. The algorithm calculates a linear interpolation 
for missing pixels with four neighbors, then three, then 
two. The order in which this interpolation is done cre- 
ates some bias in the interpolated pixels, but the general 
properties of the emission on spatial scales larger than 
the lO'-beam are not affected. 

2.2. GBT 

The goal of the GBT observation was to confirm the 
results of the HCRO observations and study the gas con- 
ditions in more detail. We used the GBT Spectrometer in 
a 7 hour session in August 2005. The GBT observations 
were more limited in spatial coverage than the HCRO 
survey; the locations are shown in Figure[5]and described 
in Table CD The locations named "GCL3" and "GCL4" 
are at the two bright peaks in the HCRO map of H109<x 
These regions were mapped in a 3x3, Nyquist-sampled 
(1' near 5 GHz) grid, with on-source integration times of 
90 s per pointing. Diagonal strips were observed starting 
at these positions and going to the Galactic northeast to 
study gradients in the velocity field. A sparse horizontal 
strip was observed at b = 0?45, to sample the changes 
in the east-west direction. Finally, two positions outside 
the GCL, "GCL1" and "GCL7", were observed to mea- 
sure background line emission. The strips and individual 
pointings had integration times of 60 s per pointing. 



The spectrometer was configured with four windows 
of width 200 MHz observing eight Ha and Hea tran- 
sitions from n =106 to 113. Each pointing produced 
dual-circular polarizations of four 8196 channel spectra. 
The spectral windows, centered near 5.37, 5.08, 4.81, and 
4.56 GHz, also covered other transitions, including: H/3 
(7 transitions with n =134-142), Hy (8; n =152-162), 
RS (8; n =167-178), and He (8; n =180-191). The spec- 
tra had 24.4 kHz channels, equivalent to 1.5 km s _1 at 
4.9 GHz. All velocities in this work are given in the 
local standard of rest (LSR). Calibration for each scan 
was done by position switching to a position two degrees 
north in Galactic latitude. 

Calibration and analysis was done using tools released 
with GBTIDL vl.2.1 and custom IDL programs. Spectra 
were calibrated using the standard (on-off)/off method. 
The typical line analysis would apply standard calibra- 
tion, average both polarizations and lines of the same 
transition (e.g., a ll H/3 transitions; rest frequencies from 
iLilley fc Palmerl |l968)), subtract a third order polyno- 
mial baseline fit, and fit a Gaussian line profile. It is 
worth noting that the standard error propagation used 
does not account for the uncertainty in fitting a baseline 
to the spectra, and are likely to underestimate the true 
errors slightly. By fitting with various assumptions, we 
estimated the baseline-induced uncertainty in the line fit 
is always less than the standard statistical errors in the 
line fit and generally decreases with signal to noise ratio 
of the line. 

All spectra presented here are corrected for the beam 
efficiency and atmospheric opacity, and thus are in units 
of main beam brightness temperature. The beam effi- 
ciency is assumed to be equal to 89%. 6 This work 
corrects for mean atmospheric absorption opacity of 
8.6xl0 -3 , giving a correction factor of 2-5% for these 
observations. 

The GC region is seen at low elevations from Green 
Bank, which makes the calibration of the continuum 
levels difficult. Unfortunately, the standard continuum 
calibration method could not properly account for at- 
mospheric emissi on. Since the observations described in 
lLaw et ail (|2008f ) were made with a calibration method 
optimized for the study of continuum levels, we use re- 
sults from that observation when continuum fluxes are 
needed. 

3. RESULTS 
3.1. HCRO 

Figure Q] shows the HCRO line intensity and velocity 
compared to 5 GHz continuum emission. The radio re- 
combination line is brightest along two vertical ridges 
near (Z,&)=(0?0, 0?5) and (359?4, 0?5). There is also 
a slight increase in brightness at the top of the line- 
emitting region, near (Z, 6)=(359?6, I?0), like a cap to 
the structure. The longitude of the line emission north 
of the plane is centered near I = 359?6 and the central 
longitude shifts slightly to the west at higher latitudes. 
The overall morphology of bright limbs with a cap mir- 
rors that of the radio continuum emission. 

The right side of Figure Q] shows the best-fit velocity 
of the H109a line. Two striking characteristics of the 
line velocities in this region are their small values and 

6 See http : //wwwlocal . gb . nrao . edu/gbtprops/man/GBTpg/GBTpg_tf . html 
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the absence of any simple pattern or gradient. North 
of b = 0?3 the line velocity does not exceed the range 
of -5 to +5 km s _1 . This is consistent with previ- 
ous o bservations of radio recom b ination lines in the re- 
gion (lLockman fc Gordonl Il973t iPauls fc Mezgeri 119751 : 
Anis h Roshi fc Ananth aramaiah 1993) ■ 

There is significant emission at every HCRO pointing 
between the two ridges of emission. At the lowest, the 
line brightness ranges from 10-30 mK, while the brightest 
line emission is 72 mK, at the peak of the western ridge. 
The emission outside of the ridges of emission is generally 
less than 10 mK. Detailed analysis of the gas properties 
depends on parameters that are not well constrained by 
the HCRO data. The GBT results are more constraining 
and are consistent with the HCRO results, so the results 
are given below in § 13.2.31 

3.2. GBT 

3.2.1. Morphology 

Figure [3] shows plots of the average H106 — 113a line 
brightness for the three strips of observations. The hor- 
izontal strip at b = 0?45 is similar to the HCRO data, 
with two recombination line peaks near the radio con- 
tinuum peaks. The eastern peak of the line emission is 
offset west the continuum peak by roughly 0?15. The 
western peak of the line and continuum emission fall at 
the same position (I ~ 359.35 — 359.4), within the 5' spa- 
tial sampling of the line observations. This is consistent 
with the morphology observed by the HCRO and again 
shows how the structure is similar to the radio continuum 
emission. 

As with the HCRO observations, the line brightness 
between the two emission peaks is significantly brighter 
than outside of it. The GCL1 and GCL7 pointings sam- 
pled the emission toward lines of sight outside the ra- 
dio continuum emission of the GCL. The best-fit Gaus- 
sian to the average H106 — 113a line in the GCL1 po- 
sition had (T h v, Av) = (0.0103 ± 0.0016 K, -9.5 ± 2.2 
km s-\29.2 ± 5.2 km s" 1 ). No H106 - 113a line was 
detected toward GCL7 with ale upper limit of 5.9 mK. 
The line brightness for the weakest lines associated with 
the GCL have T x ~ 25 mK. 

The GCL3 and GCL4 diagonal strips of pointings, like 
the b = 0?45 strip, show a decrease in peak line brightness 
as they move away from the two line emission peaks. The 
longitude half-peak width of the east and west peaks are 
roughly 0?2. 

3.2.2. Velocity Structure 

At the brightest parts of the recombination line emis- 
sion, the line velocities observed by the GBT are consis- 
tent with those of the HCRO. Figure 2] shows how the 
average H106 — 113a line properties change as a func- 
tion of Galactic longitude for the strip at b — 0?45. One 
interesting trend in Figure |4] is the east-west asymme- 
try in the line velocities, with positive velocities on the 
east side and negative on the west. Averaging over the 
four easternmost and westernmost scans in the b = 0?45 
strip gives spectra with a single Gaussian component 
with mean velocities of 2.3 ±0.2 and —2.4 ±0.2 km s , 
respectively. This velocity gradient is consistent with the 
upper limit provided by the HCRO observations. 

Comparing the peak line brightness and line widths 
in Figure [4] shows an inverse relationship: the brightest 



line emission has the narrowest widths and vice versa. 
The top panel shows the line is brightest near I = 0° 
and I = — 0?7; the bottom panel shows that the line 
is narrowest at these positions. Since the line width is 
a useful constraint to the gas properties, we needed to 
understand the origin of this correlation. 

Two potential explanations for the correlation are that 
(1) the line width is intrinsically tied to its brightness, 
or (2) the line has two components that can affect the 
best-fit line width according to their relative strengths. 
Figure [5] shows the average H106 — 113a profile for four 
scans located in the center of the GCL that, separately, 
have best-fit line widths of 20-30 km s -1 . The aver- 
age spectrum shows a slightly more complex line pro- 
file than seen near the brightest recombination line emis- 
sion. If the profile is fit with two Gaussians, their best-fit 
components are (Tj, v, Av) = (22 ± 3 mK, 5.2 ± 0.5 km 
s" 1 , 9.3 ± 1.5 km s^ 1 ) and (25 ± 2 mK, -5.0 ± 1.3 km 
s _1 , 30.6±1.9 km s _1 ). The narrow component is similar 
to that observed near the bright ridges of line emission, 
but its amplitude is comparable to the wide component, 
such that a single-Gaussian fit is moderately wide. In 
averages of other scans, there is a tendency for a simi- 
lar, wide (20-30 km s" 1 ), low-level (10-20 mK) line to 
appear. The amplitude of this line is larger than the 
~ 10 mK seen in the background outside the GCL, so 
the wide line seems to be associated with the GCL. Re- 
gardless of its origin, it seems that the changes in best-fit, 
single-Gaussian line width shown in Figure |4] are likely 
due to the relative strength of the narrow line and this 
wider component. Thus, there seems to be narrow (9 — 14 
km s _1 ) recombination line emission for all lines of sight 
through the GCL. 

Figure [6] shows an alternative way of visualizing the re- 
combination line velocity structure using l-v diagrams. 
Three strips of observations show the change in the line 
velocity as a function of Galactic longitude. These plots 
show clearly how the brightest emission is always within 
20 km s _1 the rest velocity, and that the lines are typ- 
ically about 10 km s _1 wide. The plots also help show 
how the line profile has multiple components inside the 
GCL, for I = -0?2 to -0?45. In particular, the GCL4 
l-v diagram shows two diagonal structures near I = 0?4, 
which indicate that there are two narrow components to 
the line profile with velocities approaching ±10 km s _1 
toward the eastern end of the strip. 

To compare the emission at the edges to the center, 
an average of the H106 — 113a spectra for scans at the 
edges of the GCL3 and GCL4 strips was made. Figure 
[7] shows how the average of three scans from the eastern 
edge of GCL4 confirms the visual impression from the 
l-v diagram: the best-fit spectrum has two distinct com- 
ponents with (T h v,Av) = (13.8 ±1.4 mK, 13.0 ±1.6 km 
s" 1 , 12.5 ± 1.5 km s" 1 ) and (38.5 ± 1.6 mK, -6.4 ± 0.2 
km s _1 ,8.8 ± 0.4 km s -1 ). A similar average over the 
four easternmost scans of the GCL3 strip (with I — 
0?083 — 0?133) shows two Gaussian components to the 
average profile. The two components have (Ti,v, Av) = 
(26.9 ± 1.0 mK, 7.6 ± 0.3 km s" 1 , 17.2 ± 0.8 km s" 1 ) and 
(11.1 ± 1.4 mK, -18.3 ± 0.5 km s-\8.6 ± 1.3 km s" 1 ). 
Both components are narrow and brighter than the back- 
ground, so the emission is likely to be associated with the 
GCL. 

Similar to the HCRO observations, the GBT observa- 
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tions find no emission with velocities |u| > 20 km s _1 . 
A deep spectrum was made by averaging observations at 
the peak GCL line emission, called GCL3 and GCL4 in 
TableQ] the total on-source integration time for this deep 
spectrum was 27 minutes. The rms noise in the baseline 
level of the deep GBT spectrum gives an upper limit to 
the peak line brightness of about 1 mK for v = ±500 
km s _1 and 2 mK for v = ±1500 km s _1 The peak line 
brightness limits are 1% and 2% of the peak at v ~ km 
s _1 , respectively. 



3.2.3. Line Widths and Line Ratios 

To better constrain the gas conditions, we studied sev- 
eral recombination lines in the deep integration toward 
the peak GCL line emission. The relatively small ve- 
locities observed in these integrations makes it possi- 
ble to average them together. The two 3x3, Nyquist- 
sampled patterns near the brightest regions of the GCL 
were found to have not only similar line velocities and 
line widths, but also similar line ratios, suggesting that 
these regions have similar physical conditions. 

Table [2] lists the line measurements and Figure [8] shows 
the highest and lowest significance fits in the deep spec- 
trum toward the GCL. A Gaussian line is fit to the av- 
erage line profile for all transitions of a given An (e.g., 
An = 1 for "a"), which includes 7 or 8 transitions in 
the range given in Tabic [2j Generally, these averages 
over several n values are treated like a single transition 
with n and v w 4.95 GHz, or the mean of the transitions 
averaged. 

One of the most surprising characteristics of the lines 
shown in Tabic [2] is their unusually small widths. The 
average H106 - 113a, H134 - 142/3, H152 - 162 7 , and 
H167 — 178(5 transitions have widths that are simi- 
lar to each other and have an error-weighted mean of 
Av = 13.5 ± 0.2 km s _1 . The average H180 - 191e 
line is significantly wider and is discussed in more de- 
tail below. Typical H II regions ha ve much larger line 
widths, with Ay 20 km s" 1 (jShaver et all 119831 : 
lAfflerbach et al.lll996f h A hydrogen recombination line 
width of 13.5 ± 0.2 km s _1 i s equivalent to a Dopple r 
temperature of 3960 ± 120 K (|Roelfsema fe Gosslll99l . 
which is also significantly less than values seen in typical 
H II regions (T D « 5000 - 10000K). The Doppler tem- 
perature represents the combined effects of gas motion 
and thermal broadening, so it is an upper limit on the 
electron temperature. The line width is roughly constant 
for transitions from H106 — 113a to H167 — 1785, which 
represent transitions with mean n ranging from 109.5 to 
172.5. 

The narrow line widths with central velocities 
near km s _1 are highly suggestive of stimu- 
lated emission by ioniz ed gas along the line of 
sight to the GC region (ICersqsimo fc Magnanil 119901 : 
lAnish Roshi fe Anantharamaiahlll997f ). The present ob- 
servations can test for the possibility of stimulated line 
emission by comparing the line ratios of transitions with 
different n. All transitions are observed with a simi- 
lar beam size, which means that they are probing the 
same volume of space and c an be meaningfully compared 
(Du pree fe Goldberdfl970h . In this case, the integrated 
line ratio between transitions to states n and m takes on 



the form 



R » i?LTE 



T e T c f3 n 



(1) 



where -Rlte = n 2 /n,n+An/w 2 /m,m+Am, b and f3 are 
the departure coefficients for non-LTE effects, and 
T e and T c are the electron and background contin- 
uum brightness temperatures (iDupree fc Goldberelll970l ; 
iBrocklehurst fc Seatonl Il972t iRoelfsema fc Gossi I1992T . 
The departure coefficients vary with n e , T e , and n, 
and have been calculated for many t ransitions un- 
der a wide range of gas conditions (jDupred Il972t 
ISalem fc Brocklehurstlll979l ). To estimate non-LTE ef- 
fects, we calculated departure coefficients for T e = 
20 - 10000 K and n e = 1 - 1000 cm" 3 . The lowest 
temperature and density used in the calculations cor- 
respond to the conditions that may cause stimulated 
emission observed toward the GC region at 1.4 GHz 
(jLockman fc Gordonlll973t) . The LTE line ratios depend 
on the oscillator strengths for the transit ions, which ha ve 
relative values known to high precision (jDupr cc 1969). 

Table [3] shows the integrated line ratios and the pre- 
dictions for LTE and non-LTE conditions. The ratios 
show excellent agreement with LTE for all transitions 
observed. The predicted line ratios for other models do 
not agree with the observed line ratios. 

One model, with T e = 1000 K and n e = 1000 cm" 3 , 
predicts line ratios similar to the observed and LTE val- 
ues. The amount of stimulated emission can be calcu- 
lated as 

(2) 



r i /r ; LTE «6„[i-/?„(T c /r e )], 



where TJ is the line brightness, and th e background con - 
tinuum brightness is T c w 1 - 2 K (|Law et al.ll2008h . 
For T e = 1000 K and n e = 1000 cm" 3 , the amplifica- 
tion of the line is only a few percent. A much stronger 
background continuum is required to cause significant 
stimulated emission, with 10% amplification for T c ~ 25 
K; this is more than 10 times brighter than the observed 
continuum in the region. Thus, the line ratios severely 
constrain the contribution of stimulated emission to the 
recombination line emission for the brightest emission 
observed toward the GCL. 

There are two other points that suggest that the recom- 
bination line emission is not stimulated. First, if back- 
ground radio continuum was stimulating the recombina- 
tion line emission, the morphology of the line emission 
would closely follow the background continuum emissio n 
(|Lockman fc Gordonlll973l : ICersosimo fc Magnanil[l9 90). 
The HCRO and GBT observations show that the peak 
line emission, particularly in the eastern half of the GCL, 
is significantly offset from the peak continuum emission. 
The similar line width, line velocity, and line ratios to- 
ward the eastern and western recombination line peaks 
suggest that the entire structure emits by a similar mech- 
anism. Second, the width of the average H180— 191e line 
is significantly larger than that of the other lines, which 
is inconsistent with the idea that it is stimulated. How- 
ever, the H189 — 191e line represents the highest n states, 
where stimulated emission is expected to be strongest. 
These points strengthen the case for an LTE origin for 
the recombination line emission, which implies that the 
narrow line widths are related to intrinsic properties of 
the ionized gas. 
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Using the LTE assumption, the line properties give use- 
ful constraints on the intrinsic gas properties. As men- 
tioned before, the width of the average H180 — 191e line is 
significantly larger than the other transitions. Widths of 
lines that are not stimulated are generally dominated by 
thermal and Doppler broadening, although these effects 
are not expected to change between the different tran- 
sitions observed here. The most common n-dependent 
broadening effect is collisional broadening (a.k.a. "im- 
pact" broadening), in which inelastic e lectron collisions 
preferentially broaden lines wit h large n (|Lang fc Willsonl 
ll97llRoelfsema fe Goss|[l"992T) . Expressing the total line 
width as the quadrature sum of all broadening effects 
gives: 



km s 



4.31 / n y 
~A~~~ V100/ 



10 4 cm- 3 V 10 4 K 



T 



(3) 

where Avd is the Doppler line width. 7 

The strong dependence of collisional broadening on n 
can be used to constrain n e , since Av^ is known and the 
dependence on T e is weak. Equation [3] was fit to the 
distribution of Hydrogen line widths assuming Avd = 
13.5 km s , which is the error- weighted mean of the 
H106 - 113a to H167- 1785 line widths 8 and T e = T d = 
3960 K. Figure [9] shows the best-fit line width model with 
n e = 950+310 cm~ 3 (ler). The fit quality is relatively 
poor, with x 2 l v = 14.4/4, which suggests that the errors 
may be underestimated somewhat. For a more realistic 
fit, we model additional uncertainty in the line widths as 
equal to the statistical uncertainty for the H180 — 191e 
line and scaled by the line peak signal to noise ratio; 
this uncertainty is consistent with variation in the fit line 
width with changes in the baseline fit. Adding this line 
width uncertainty in quadrature to the statistical error, 



910 



-310 
-450 



with 



the best-fit electron density is n e 

a more realistic xl = 8.4/4. 

The Hel06 — 113a line was detected in the deep spec- 
trum with 18cr significance. The ionized helium abun- 
dance is measured by the integrated line ratios Y + = 
Inca/Ina- In the GCL, this ratio is 0.075 ±0.007, simila r 
to the solar value (|Landll974t iRoelfsema fe Gosslll992l) . 

Typically, the line-to-continuum ratio for a radio re- 
combination line is a useful way to measure the electron 
temperature of the emitting region. The LTE electron 
temperature is: 



T, 



LTE 



6943 



11 T r 



1 



GHz T t Avl + Y+ 



0.87 



K. 



(4) 



which is related to the actual electron temperature by: 

n -0.87 



nLTE 



1 



(3 n T c 



(5) 



where t c is the optical depth of the background contin- 
uum emission. However, much of the continuum emission 

7 Technically, collisional broadening takes the shape of a Lorentz 
distribution, which, when convolved w ith a Gaussian, produces a 
Voigt profile (Roelfscma & Goss 1992). However, at the low sig- 
nal to noise ratio observed for the H180 — 191e transition, no line 
wings are apparent and a Gaussian is a close approximation to the 
expected line shape. 

8 Allowing Adj to vary did not significantly change the fit values 
or its quality. 



in the GCL is nonthermal (|Law et al.ll2008D . so the line- 
to-continuum ratio is an upper limit on the electron tem- 
perature. For the deep spectrum, the continuu m flux was 
estima ted from the 5 GHz map described in lLaw et aTl 
;'2()(),S; to be r c w 0.8 K, which gives TiAv/T c = 2.0 
km s _1 , for the averaged H106 — 113a line. This line- 
to-continuum ratio gives to an upper limit on the LTE 
electron temperature of T^™ < 5220 K, which is con- 
sistent with the upper limit on T e from the narrow line 
widths. For the largest non-LTE effects that are consis- 
tent with the observed line ratios (n e = 100 cm~ 3 and 
T e = 1000 K), T e = (0.855) a87 T e LTE < 4550 K, but for 
higher densities and temperatures, T e is within a few per- 
cent of T e LTE . Alternatively, we can use the observed line 
width to constrain the electron temperature. The upper 
limit to the electron temperature gives a lower limit on 
the line-to-continuum ratio of TiAv/T c > 2.75 km s _1 
and an upper limit on the thermal continuum of 0.6 K. 
Thus, at most about 70% of the radio continuum toward 
the recombination line peaks of the GCL is thermal. 

The upper limit to the electron temperature constrains 
the emission measure. Assuming r l <C 1 and LTE for the 
average H106 — 113a transition, the continuum emission 
measure is: 

EM C = 8.5 x 10" 3 (1 + Y+)AvT 1 T 3/2 e- x "pc cm" 6 , (6) 

where At; is in units of km s -1 and e~ Xn f=a 1, for 
T e > 1000 K (jDupree fe Goldberd[l970T h Using the av- 
erage H106 — 113a line parameters for the deep spec- 
trum gives EM C « 3850(T e /3960 K) 3 / 2 pc cm" 6 . The 
eastern and western peak line brightnesses are some- 
what different, giving EM C « 3080(T e /3960 K) 3 / 2 and 
4570(T e /3960 K) 3 / 2 pc cm" 6 for the east and west, re- 
spectively. 

4. DISCUSSION 
4.1. The Structure of the Ionized Gas 

The GBT and HCRO observations have found radio 
recombination line emission throughout the radio con- 
tinuum emission of the GCL. The strongest emission is 
found along an eastern and western ridge, but significant 
emission is found between these ridges. All of this emis- 
sion has an unusually narrow line width and LTE line 
ratios, suggesting that it is all a part of the same struc- 
ture and distinct from background Galactic emission. 

The morphology of the ionized gas suggests that it 
is associated with the radio continuum emission of the 
GCL. The two ridges and cap of recombination line emis- 
sion are parallel and adjacent to simi lar structures in 
radio continuum maps (|Law et al.ll2008h . If the radio re- 
combination line and continuum emission are colocated, 
the line emission may be in the G C region, as has been 
argued for radio continuum GC L (jLasenbv et alj Il989l : 
iBland-Hawthorn fc Cohenl 120031 ) . A more detailed ar- 
gument is left for a future paper; for the results shown 
below, we derive physical parameters for the ionized gas 
assuming a distance of 8 kpc. 

The morphology of the line emission is consistent with 
a three-dimensional shell model. In this model, the two 
ridges of emission ar e the shell edges, where the col - 
umn density is largest. [Bland-Hawth orn fc Cohenl (|2003l ) 
model the MSX mid-IR emission from this region as a 
"telescope dome" shell specified by a height, radius, and 
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shell thickness. Assuming this geometry, the radius of 
the dome is roughly 0?3 and the height is roughly 1°. The 
GBT observations measured the ridge half-peak thick- 
ness to be about 0?2; simple simulations suggest that for 
a shell geometry, the true shell width is about half the 
apparent thickness, or 0?1. Assuming this shell is in the 
GC region, the dome radius is 40 pc, the total height is 
140 pc, and the width is 15 pc. The path length through 
the shell edge would be about 50 pc. Combining this 
with the EM measurement, gives an rms electron den- 
sity, VT^I) ~ 8.8(T e /3960 K) 3 / 4 cm" 3 . 

The morpho logy of the G CL is reminiscent of "Galac- 
tic worms" (jKoo et al.l ll992). Worms are ionized cavities 
formed when supernovae and stellar winds from H II re - 
gions blow out of the Galactic disk (Hei les et al.lll996f) . 
A difference between the GCL and a worm is that the 
GCL has a cap at its top. The GCL may be thought 
of as a frustrated worm, an ionized cavity that has not 
yet blown out of the Galactic disk. A degree-scale struc- 
ture like this woul d not have been detected by previous 
surveys for worms (jKoo et al.lll992h . 

4.2. Properties of the Ionized Gas 

It is critical to show that the recombination line emis- 
sion is in LTE if we wish to derive intrinsic gas conditions. 
Previous observations in the GC region found narrow 
recombination lines, particularly at lower frequencies, 
and co ncluded that the emission was most likely stim- 
ulated (lLockman fc Gordonlll973|;lPauls fe Mezgen ll975l; 



lAnish Roshi fe Anantharamaiahl 1997f ). However, there 
are several characteristics of the line emission observed 
here that show that the gas is in LTE, including well- 
constrained line ratios consistent with LTE, the lack of 
coincidence between background continuum brightness 
and line brightness, and evidence for collisional broaden- 
ing in the line widths. 

Assuming the emission is not stimulated, the Doppler 
widths of these lines provide a strict upper limit to the 
electron temperature of 3960 ± 120 K. This limit to the 
electron temperature is among the lowest values observed 
in the Galaxy, although radio recombination line obser- 
vations of other H II regions have found electron and 
Doppler temperatures as low as 4000 K (|Shaver et al.1 
1198a lAfflerbach et al.lll996l M. Goss 2006, private com- 
munication). Some of the best-fit line widths discussed 
in this work are as low as 9 km s , which corresponds to 
a Doppler temperature of ^1800 K, much smaller than 
observed elsewhere in the Galaxy. 

The dominant effect in determining the electron tem- 
perature of ionized gas is the cooling efficiency; high 
metal abunda nce allows ionized gas to cool quickly 
(jMathisI 119851 ). Studies of recombination line emission 
from H II regions have used this effect to establish that 
there is a metallicity gradient with distance from the GC 
(jShaver et all Il983t lAfflerbach et al.lll996h . The gradi- 
ent is seen as a trend in T e , with coolest and most metal- 
rich gas in the GC re gion. The relation observed by 
lAfflerbach et~aTl (Il996h pre dicts T c 5500 K in the GC 
region, while IShaver et all (|1983l ) predicts T e m 3100 K. 
Regardless of the exact value, these relations suggest that 
the ionized gas in the GCL is located near the GC, sup- 
porting the morphological connection. 

The broadening of the line width with increasing elec- 
tronic state, n, is best explained by collisional broad- 



ening. The best-fit model to the line widths gives 
n e — 910^450 cm" 3 . The best- fit density is consis- 
tent with the lack of stimulated emission, which would 
change the line ratios for n e < 100 cm -3 . Using the 
rms and true electron densities, we calculate the volume 
filling factor for the deep spectrum as / = (n 2 )/n 2 = 
(9±f ) x 10" 5 (T e /3960 K) 13 (50 pc/L) cm" 3 , where L is 
the path length through the shell, assuming a distance 
of 8 kpc. Since the emission measure scales as (rig), the 
derived filling factor applies to the densest gas along the 
line of sight. This filling factor is much smaller than 
that seen in the Galactic disk (/ ~ 10~ 2 ; Heil es et al.l 
fl99l . but compar able to that seen in G alactic outflows 
(/ ~ 10" 3 - !Q- 4 : lHeckman et al.lfl990h . 

The thermal pressure implied by the constraint on T e 
and measurement of n e is P/k = 7.2 x 10 6 (T e /3960 K) K 
cm" 3 . This pressure is about 10 times larger t han the 
total gas pressure near the Sun (lBloemenlll987f). but it 
is not unusua l for t h e GC region (ISpergel fc Bln9 [19921: 
iMartin et all l2004t iKovama et all 119961 : iMuno et al l 
l2004h . This suggests that the ionized gas is in the GC 
region and in equilibrium with its molecular and hot, 
x-ray-emitting gas. 

4.3. Mass and Ionization of the Gas 

Assuming a shell geometry and GC distance to the 
recombination line emission allows us to estimate other 
physical properties of the gas. For ionized gas in LTE, 
the mass is parameterized as: 

M = 0.419(T e /10 4 K) ' 175 (S 5 GHz/Jy) ' 5 ( J D/kpc) 2 - 5 (e G /arcmm) 

(7) 

where D is the distance to the gas, 9q is the an- 
gular size of the emission, and the geometry is as- 
sume d to have a cylidrical sh ape (|Pauls fc Mezgerl 
[1971 iPanaeia fc W almslcvi fl978h . Since the gas is 
in LTE, the line-to-continuum ratio is T/Aw/T c = 
2.75(T e /3960 K)" 87 km s" 1 , as shown in § l3~Q The 
integrated line intensity over the ionized shell of gas in 
the HCRO observations is about 1.2 K* 13.5 km s = 16 
K km s _1 , for b > 0?2 and subtracting 10% for the back- 
ground. This gives a 5 GHz continuum flux of 50 Jy and 
a ionized mass M = 2 x 10 5 (T e /3960 K) 61 M Q in the 
shell. Uncertainties in the mass estimate are probably 
dominated by incomplete sampling of the HCRO survey 
and are about 20%. 

The mean electron density constrains the number of 
Lyman continuum photons required to ionize the gas. 
The Lyman continuum photon flux is N]_, y = Vn e naa^ , 
where V is the volume of the emitting region and 
a< 2 > » 6.6 x IP" 13 cm 3 s" 1 , for T e = 3 x 10 3 K 
(Schcfficr & Elsasserl ll988l ). Using the integrated ther- 
mal continuum flux density, we find that the gas requires 
an ionizing flux N Ly = 7 x 10 49 (T e /3960 K) 122 ) s" 1 . 

The estimated ionizing flux is about seven times the 
Lym an continuum flu x of an 07 V star of about 10 49 
s -1 (|Smith et al.ll2002f ). Assuming it is a shell in the GC 
region, the source of ionizing photons would most likely 
be stars in the plane. Based on geometry, at most half 
of the stellar flux would reach the shell, so the required 
ionizing flux is equivalent to more than 14 07V-type stars, 
or 1.4 x 10 50 s" 1 . The Arches, Quintuplet, and Central 
star clusters, all located inside the GCL in projection, 
have measured Lyman continuum fluxes of 10 51 , 10 50 9 , 
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and 10 50 - 5 s _1 (jFiger et al.lll999f) . Even if extinction is 
large, these clusters could easily ionize the gas associated 
with the GCL. 

4.4. Kinematics 

The observations confirm previous observations that 
found line velocities ne ar zero toward the GCL 
(|Lockman fc Gordon! [l973h . However, the GBT l-v di- 
agrams show that the line velocities tend to be much 
more complex, particularly between the bright ridges of 
emission, from I = 0?0 to — 0?5. Unusually narrow lines 
are brightest near the GCL and have physical conditions 
most consistent with the GC region. This narrow line 
observed across the GCL at b = 0?45 shows a velocity 
gradient similar to Galactic rotation, indicating that the 
gas may be connected to the disk. 

The low velocity of the ionized gas has been used 
to argue that it is in the foreground to the GC region 
(|Pauls fc Mezgerll975f) . However, recent work has shown 
that a significant fraction of gas in th e central few hun - 
dred parsecs has velocities near zero. lOka et all {2005) 
studied H% absorption lines toward the GC region; these 
lines are excited at temperatures and densities found in 
the central few hundred parsecs and have only been ob- 
served toward the GC region. They found that roughly 
1/3 of column density has velocities near zero. Thus, 
lo w velocity ga s is quite common in the GC region. 

Soful (|1996f ) argued that the GCL is associated with 
molecular gas rotating about the GC at ± ~ 100 km 
s _1 . While that association relies on a positional co- 
incidence in a complex region, it is worth considering 
how it would relate to the ionized gas. In particular, 
how can the large difference in velocity of these two 
components be explained? The motion of this molec- 
ular gas is best explained by the gas dynamics in a 
barred gravitational potential (jBinnev et al.lll99lh . Un- 
der this model, the gas possibly associated with the 
GCL is either in or transitioning to th e inner "x2" or- 
bits (|Contopoulos fc Mertzanideslll977T) . In light of this, 
the low velocity of the ionized gas in the GCL may be 
tied to the shocked molecular gas on these transitioning 
orbits. 

Alternatively, the low velocity of the ionized gas may 
indicate that it as been decelerated by an ambient mag- 
netic field. Equating gas ram pressure to a magnetic 
pressure allows us to estimate the magnetic field re- 
quired to stop the gas on a dynamic time scale . For 
the m olecular gas associated with the GCL by ISofud 
(1996), these pressures are equal for B = 1 mG. The 
possibility of large-scale, mG-strength fields in the GC 
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region is being active ly debated ([Yusef-Zadeh et al]|1987t 
iLaRosa et aUl2005t iBoldvrev fc Yusef-Zadehl l2006h : a 
weaker field would decelerate more slowly, but could ul- 
timately have the same effect. 

5. CONCLUSIONS 

We have presented analysis of new and archival radio 
recombination line observations toward the GCL, a radio 
continuum shell believed to be evidence of a mass outflow 
from our GC region. The observations have found that 
the radio line emission has a morphology strikingly sim- 
ilar to the radio continuum, strongly arguing that they 
are associated. The line and continuum emission appears 
as a limb-brightened shell. A simple shell model for the 
structure is best fit with a radius of 0?3, a height 1°, and 
a shell width 0?1. 

Recombination lines are detected in averaged profiles 
covering transitions from H106a up to H191e. Diagnos- 
tics derived from the detected lines show that the emis- 
sion is not stimulated, allowing us to constrain the elec- 
tron temperature and gas density from the line widths. 
The electron temperature is unusually low, with an up- 
per limit of roughly 4000 K. This temperature limit, and 
the gas pressure, are consistent with conditions in the 
GC region. 

Assuming that the gas is organized as a shell and 
is in the GC region, we derive a mass of M = 2 x 
10 5 (T e /3960 K) 61 M Q . This mass is similar to that ob - 
served in dwarf starburst outflows (|Veilleux et al.| [2005) . 
The filling factor of the ionized gas is also similar to that 
seen in starburst outflows. In the canonical model for 
such outflows, the emission comes from the ionized sur- 
face of clumps of mole cular gas entrained in the outflow 
(jHeckman et al.lll990t l. 

Modest observations with radio interferometers cur- 
rently in development could greatly expand on the results 
presented here. A simple estimate using the filling fac- 
tor and the size of our deep integration shows that these 
chimps could be resolved at size scales of ~9". Deeper 
observations of more transitions would strengthen the 
significance of collisional broadening and detect spatial 
structure in the electron density and gas pressure. Fi- 
nally, high-spectral resolution observations would map 
the line velocities and help us understand its dynamics. 
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TABLE 1 

GBT Radio Recombination Line Pointings 



Name 


(dag) 


(deg) 


Map type 


Exposure 

00 


Spacing 
(arcmin) 


GCL3 


0.0 


0.5 


3x3 grid 


90 


1 


GCL4 


359.4 


0.5 


3x3 grid 


90 


1 


GCL3 strip 


+0.016-0.133 


0.516-0.633 


diagonal strip 


60 


1.4 


GCL4 strip 


359.416-359.666 a 


0.516-0.766 


diagonal strip 


60 


1.4 


b = 0.45 strip 


0.25-359.25 


0.45 


horizontal strip 


60 


5.4 


GCL1 


0.3 


0.5 


single point 


60 




GCL7 


359.1 


0.8 


single point 


60 





One pointing, near I = — 0?48, was accidentally skipped. 



TABLE 2 

Detected Lines and their Properties in Deep Integration Toward GCL 
Recombination Line Peaks 



Transition 


Ti 


"LSR 


A.v 


jobs a 


oobsb 




(K) 


(km s" 1 ) 


(km s — 1 ) 


(K km s" 1 ) 





H106 - 113a 
Hel06 - 113a 
H134 - 142/3 
H152 - 1627 
H167 - 1785 
H180 - 191e 



0.1257 ±0.0024 0.71 ±0.13 13.45 ± 0.30 1.802 ±0.052 1 

0.0110 ± 0.0006 0.97 ±0.32 11.50 ±0.74 0.135 ±0.012 0.075 ± 0.007 

0.0339 ± 0.0008 0.75 ±0.15 13.72 ± 0.36 0.495 ± 0.017 0.275 ± 0.012 

0.0144 ± 0.0007 0.36 ±0.35 15.26 ±0.83 0.233 ± 0.017 0.130 ±0.010 

0.0091 ± 0.0006 1.98 ±0.40 12.47 ± 1.05 0.120 ±0.013 0.067 ± 0.007 

0.0030 ± 0.0003 1.67 ± 1.34 24.47 ±3.15 0.073 ± 0.012 0.041 ± 0.007 



/r,n 



1.065TiA« b ii° ba = I° b °/I° 



TABLE 3 

Comparing Line Ratios to LTE and Stimulated Emission Models 



Transition 



R 20, 
H a 



D 100,10 
H a 



D 100,100 
H a 



D 100,1000 
H a 



D 1000,100 
H a 



D 1000,1000 
H a 



□10000,10 
H a 



R 10000,100 R 10000,1000 
Ha Ha 



H106 - 113a 1 

H134 - 142/3 0.275±0.012 

H152 - I627 0.130±0.010 

H167 - 1785 0.067±0.007 

H180 - 191e 0.041±0.007 



1 

0.279 
0.127 
0.073 
0.047 



1 

0.322 
0.224 
0.146 
0.112 



1 

0.381 
0.198 
0.118 
0.077 



1 

0.300 
0.137 
0.083 
0.055 



1 

0.295 
0.133 
0.077 
0.049 



1 

0.281 
0.123 
0.071 
0.053 



1 

0.275 
0.124 
0.072 
0.048 



1 

0.417 
0.193 
0.107 
0.055 



1 

0.242 
0.098 
0.053 
0.051 



1 

0.245 
0.103 
0.058 
0.047 



Note. — The line ratios for LTE conditions correspond to the "mean" transition of the range of n values shown. Oscillator strength s were taken at 
the mean n or the average of the nearest two n values, for fractional n. Stimulated emission models labeled with (T e , n e ) in K and cm " 3 . IDupreel i 1972ft 
and Salem &. Brocklchurst (1979) were used for T < 100 K and > 100 K, respectively. 
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Fig. 1. — Left: The gray scale shows GBT 5 GHz radio continuum emission toward the GCL with contours of H109a brightness from 
the HCRO at T a =15, 20, 30, and 40 mK. Right: Same as the left panel, but with a contour of H109a line velocity at km s . The areas 
to the far west and east have negative velocities and between the lines are generally positive velocities. 



0.900 



0.800 



GCL3 strip 



GCL4 strip 




GCL7 



359,200 



Fig. 2. — Posit ions of the GBT observations toward the GCL are shown as circles on a GBT 5 GHz radio continuum survey of the region 
IILaw et al.H200a) . The circles have a diameter of 2'.5, which is the FWHM of the GBT beam at 5 GHz. 
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0.2 0.0 -0.2 -0.4 -0.6 -0.8 

Galactic Longitude (deg) 



Fig. 3. — Peak brightness temperature of Ho as a function of Galactic longitude for three strips of GBT observations through the GCL. 
The b = 0.45 strip is shown with black crosses, diagonal "GCL 4 strip" shown with red stars, and diagonal "GCL 3 strip" shown with 
blue diamon ds (see Ta ble [T] for strip specifications). Error bars show 1 a uncertainties. The line shows the GBT 5 GHz continuum flux 
described in Law ct al. (2008), which defines the extent of the radio continuum GCL. 
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Fig. 4. — Plots of the properties of the average Ha recombination line for the series of pointings across the GCL at b = 0?45. Line 
properties found by fitting a single Gaussian to the profile. The top plot shows the peak line brightness temperature, the middle plot shows 
the line velocity, and the bottom plot shows the line width. The brightest recombination line emission (near I ~ 0° and — 0?6) has the 
narrowest line width. 
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Fig. 5. — Average Ho spectrum for four scans near the center of the GCL with a line showing the best-fit model described in the text. 
Fitting one Gaussian gives a moderately wide fit, but fitting two Gaussians gives a narrow and wide component of similar amplitude. This 
shows how the assumption of a single Gaussian component can lead to a inverse relationship between line width and amplitude and that 
the unusually narrow line is observed throughout the GCL. 
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Fig. 6. — l—v diagrams of the average Ha radio line brightness for the GCL3 strip (left), b = 0?45 strip (center), and GCL4 strip (right). 
The GCL3 strip starts near the eastern peak of the recombination line emission of the GCL and goes toward the northeast. The b = 0?45 
strip spans the entire east-west extent of the GCL. The GCL4 strip starts near the western peak of the line emission in the GCL and goes 
toward the northeast. Detailed descriptions of these strips of observations are given in Table Q] The right l-v dia gram docs not account 
for the missing scan at I = — 0?48, which gives errors in I of order l',4 near that longitude. 
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Fig. 7. — Plot of the average Ha spectrum for three scans inside the GCL, on the eastern edge of the GCL4 strip. The profile requires 
two narrow Gaussians for an adequate fit, with velocities of roughly +14 and —6 km s^ 1 . The unusually narrow lines suggest that this gas 
is associated with the GCL 
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Fig. 8. — Top: Average Ho spectrum for eight transitions with n =106 to 113 in the deep average spectrum of the GCL. The best-fit 
Gaussian line described in Table [2] is overlaid. The lump near v = 5.446 GHz is the smeared average of the nearest He« transitions. 
Bottom: Same as above, but for eight He transitions with n =180 to 191. The brightest line, at u = 5.407 GHz, is H108a, which is about 
4 MHz from H183e. Note that the frequencies of the spectrum shown here correspond to that of the highest-frequency transition observed, 
but that all velocities and line properties are calculated by assuming an idealized "mean" of all transitions with a given An. 
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Fig. 9. — Plot of the distribution of hydrogen line widths observed by the GBT as a function of electronic state n. The solid line shows 
the best-fit model to the widths, considering the effects of collisional broadening; the dotted lines show the ±1<t fits. The normalization of 
the curves corresponds to an electron density n e = 910^45q cm -3 , assuming T e = 3960 K. 



